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Abstract — A theory called Adaptive Disturbance Tracking 
Control (ADTC) is introduced and used to track the Tip Speed 
Ratio (TSR) of 5 MW Horizontal Axis Wind Turbine (HAWT). 
Since ADTC theory requires wind speed information, a wind 
disturbance generator model is combined with lower order 
plant model to estimate the wind speed as well as partial states 
of the wind turbine. In this paper, we present a proof of 
stability and convergence of ADTC theory with lower order 
estimator and show that the state feedback can be adaptive. 


I. Introduction 

Large wind turbines are operated in three different 
regions called Region I, Region II and Region III, and the 
available wind speed determines the region of operation. 
Region I is the startup region and wind speed is not sufficient 
to produce the power. When wind speed becomes large 
enough to produce the power but not enough to produce the 
rated power, then it operates in Region II. In Region III, wind 
speed is normally larger than the rated wind speed and 
produces the rated power. 

The amount of power that can be captured using wind 
turbine is given by [1]: 

1 , 

P = — pAC p co ( 1 ) 

Where, P is the power, A is the area of the rotor disc, C P is 
the power coefficient, p is the air density, co and is the free- 
stream wind velocity. 

The power coefficient C P is the function of tip-speed ratio 
(! X ) and the blade pitch angle (p). The Tip Speed Ratio (TSR) 
can be expressed as: 



( 2 ) 
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Figure 1 Plot of C P vs Tip - Speed Ratio for different Blade 
Pitch Angle 

Where, Q is the rotor speed, R is the rotor radius. 

The variation of power coefficient (C p ) with TSR for fixed 
sets of blade pitch angle is in figure 1 [2] . From the figure, it 
can be seen that, for a certain blade pitch angle, there exist 
an optimum TSR for which C p is maximum. In Region II, 
the wind speed is not sufficient to produce the rated power 
so the TSR must be maintained constant at its optimum 
value such that the power coefficient is maximum, which 
ultimately maximizes the power captured by wind turbine. 

From equation (2), the TSR is a function of the rotor speed 
and the wind speed. Since variation in the wind speed have 
effect on rotor speed, the rotor speed must be varied in same 
proportion to keep the TSR constant. 

The foundation of counteracting the persistence disturbance 
in a plant was developed by Johnson [3]. Balas refined this 
idea [4], and developed a new idea of Disturbance Tracking 
Control and used for Region II control of the wind turbine. 
The motivation behind the theory of Disturbance Tracking 
Control is to make the wind turbine track the wind speed, 
which ultimately tries to keep the TSR constant at some 
optimum value. 

The TSR tracking error (e) is introduced as the deviation of 
the actual TSR (X) from the operating or optimum TSR 
(Aop). 
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where (Q op ,co op ) is the desired turbine operating point 
corresponding to the desired TSR (^ 0P ). We let the Output 
Tracking Error be 


e y =C1 T -Q* co 
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(4) 


where (Q) is the tracking ratio. 

and think of (Q), the turbine speed, as a measured output of 
the turbine and (co), the wind speed, as a disturbance on the 
turbine. Then DTC uses a feedback control law that 
produces: 


e v = Q — Q* a> — >0 

y ^ *- >o ° (5) 

This (approximately) produces tracking of the desired TSR: 
e = X-l nn — >0 [5], 

op t—t 00 L J 


II. Adaptive Disturbance Tracking Control with 
Partial State Estimation and State Feedback 


In this section we introduce further modification of the 
theory presented in [5] with addition of wind speed and 
partial states estimation, and partial state feedback. 

The wind turbine is assumed to be modeled by a linear, 
time -invariant, finite-dimensional system: 


f x = Ax+Bu p +Tu d 
[ y = Cx; x(0) = x 0 


(6a) 


and the plant model with partial states of (6a) is expressed 
as: 


All matrices in Eqs. (6)-(7) have the appropriate 
compatible dimensions. Such descriptions of persistent 
disturbances were first used in [5] to describe signals of 
known form but unknown amplitude. Equation (7) can be 
rewritten in a form that is not a dynamical system, which is 
sometimes easier to use: 

=@Z ° ( 8 ) 
[z D = L(p D 

where, (p D is a vector composed of the known basis 
functions for the solution of u D =0z D , i.e., (p D are the basis 
functions which make up the known form of the disturbance, 
and L is a matrix of appropriate dimension. The method for 
tracking persistent disturbances used in this paper requires 
only the knowledge of the form of the disturbance, the 
amplitude of the disturbance does not need to be known, 
i.e.(L,0) are unknown. In this paper, we will be interested in 
rejecting step disturbances of unknown amplitude which can 
be represented in the form of Eq. (8) as (p D =l , with (L,0) 
unknown. This has been a viable model for wind 
fluctuations in our previous work. 

Now combining equation (6b) and (7) we get a new 
augmented plant model [6]: 
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In equation (9) we used the lower order plant model to 
estimate the partial state and use the partial state feedback. 

Using the augmented plant model in (9) , a state estimator 
can be designed as: 


fx —A x +B u +r u n 
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where, the plant state, x, is an N p -dimensional vector, the 
control input vector, u p , is M-dimensional, the sensor output 
vector, y , is P -dimensional. x m is m-dimensional (m<n) 
lower order plant model state vector, y m is the P m 
dimensional model output vector. A s B C T are the state, 
input, output and disturbance matrix of plant with 

appropriate dimensions. A raj B raj C raj are the state, input, 

output and disturbance matrix of the lower order plant model 
with appropriate dimensions. The disturbance input vector, 
U D , is M D -dimensional and will be thought to come from the 
Disturbance Generator: 


[u D — 0z D 

\z D = Fz d ;z d ( 0) = z 0 
where, the disturbance state, z D , is N D -dimensional. 


(7) 


( 10 ) 

The estimator equation (10) can also be broken down into 
wind disturbance estimator and state estimator as: 


K = A J m + K u + r Jz D + K < (y P - y m ) 

(ii) 

^D= F ^D+ K D(y P ~yJ 

Using the wind disturbance state estimation, the estimated 
wind speed can be expressed as: 

u D =0z D (12) 

Our control objective will be to cause the output of the 
plant, y p , to asymptotically track some linear function of 



estimated disturbances of the form given by the disturbance 
estimator. We define the estimated output tracking error 
vector as: 

e y =y p -Qu D (13) 

To achieve the desired control objective, we want 


This aligns with the TSR tracking in Region II described 
by equation (5). 

Consider the plant given by Eq. (6a) with the disturbance 
generator given by Eq. (7) and respective disturbance and 
state estimator given by (10) and (11). Our control objective 
for this system will be accomplished by an Adaptive Control 
Law of the form: 
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Define the ideal trajectories as: 

r 

w* = Aw* + Bu* 

<y.=[C 0]w. 
e = Cw. = 0 

l y 

with 


Up = G e e y + G D (p D + G x x m 

(14) 


where G e and G D are adaptive gain matrices of the 
appropriate compatible dimensions, and G x is the state 
feedback gain matrix. 

Now we specify the Adaptive Gain Laws , which will 
produce asymptotic tracking: 


G e = -e y e y y e 
Ad — —^^dYd 
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u* = S 2 z D + G* z z* ; z* = Tw* = TS { 
=> u* = (S 2 + G*TS l )z D + G\z* 

\ j 


let 

A w = w-w* 


Ax 
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where / e ,/ D are arbitrary, positive definite matrices. Our 

Adaptive Controller is specified by Eq. (14) with the above 
adaptive gain laws Eq. (15). 

III. Stability and Convergence Analysis 


A u = u-u* 

&y = y-y* 

Ae = e —e — e 

y y y y 

A z = z — z* = Tw-Tw* =TAw 


A. Problem Formulation 

Assume a linear, time invariant and finite dimensional plant 
of (6a) and disturbance generator of (7) 

We define the estimator as: 


Aw = A Aw + B Au 


A y = [C 


0] Aw 


e = Ae = CAw 

l y y 


z = L n y + L 22 z + L 23 u 
where z = Tw + e 


now assume: 

z D =H D z or H J = [Q I D ] 


Where, 



(16) 


=> z D = H d* = H D Tw+H D e 
. • z D = z D + H D e 
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define 



\u D =&D 

\e y =y~Qu D 

i=y- QQz d 
= y~ QBz d - Q9H D e 



= e y - Q0H D e 
= (<?, ~ e y ) ~ Q&H D e 
= Ae y - QQH D e 
=>e y =CAw-Q9H D e 

adaptive Control law is: 
u = Ge y +G D (p D +G z z 

Am =u — u, —u - S 2 z D 
- u ~(S 2 z D +G]z,) 

= u - S 2 L D <p D - c/z, 

= G'e t + (G: - SL d ) <p D + G\z - G\z , 



(23) 


(24) 


(25) 


= G' (CAw - QOII e) + G (TAw + e) + AGi) 

h 

= (GC + GT)Aw + (G* - GQ0H d )e + h 
Estimator error is expressed as: 
e — z-Tw 

= [L 2l y +L 22 z +L 23 u]-T[Aw+ Bu] (26) 

[C 0]w Tw+e 

= (L 2I [C 0] + L 22 T — TA ) w + L 79 c + (Z, 23 — )w 
Choose L 23 = TB where T 3TA-L 22 T = L 2l [C 0] 

( uniquely solvable when cr(^) ncr(L 22 ) = <p) 
with L 22 chosen appropriately stable 



+ tfGjT; , B, C) is ASPR <=> CB > 0 
and P x (s) minimum phase. 



But State Feedback does not change minimum phase of a 
transfer function, proof of this will be presented later; so 
P(s) minimum phase <^> P x ( 5 ) minimum phase. 

So CB > 0 and P(s) minimum phase will guarantee that 
3G* 3 (A + BG\ C, B , C) is SPR . Consequently, we have the 
following : 


Theorem 

1) CB > 0 and P(s ) minimum phase 

2) (p D bounded 

=> Adaptive Control 
[u = G e e y + G D (p D +G z z 

[ z = L 2X y + L 22 z + L 23 u 

with gain adaption laws 


G e =-e y e T y y e -,r e >0 
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produces 
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e y = CAw-Q9H D e—^0 
e y = CAw — > 0 
with bounded adaptive gains 


This does allow adaptive G z . 


These Matrix Matching Equations can be uniquely solved 
for (iSj , S 2 ) when CB is nonsingular (which it is when 
CB>0). 

New Adaptive Controller: 


ju = G e e y + G z z + G D <p D 
| Z — L 2x C y + L 22 Z + L 22 U 


(32) 


with 

G e =-e y e T y e -,y e >0 

„ 7- (32) 
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G z =-e y z T y z \y z >0 
Use 

e y = y-Q9H D z=>z 
= L 2X e y L^z L 23 u 
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Then 

e = z — Tw = [L 2 j y 4 - L 22 z + L 23 u\ — T\Aw + Bu ] 

[C 0 ]w 7W+e 

= (L 21 [C 0] + L 22 T — TA)w + L 22 e = (L 23 — TB)u ^ 4 ^ 

Choose L 23 = TB where T 3 TA-L 22 T - L lx [C 0] 

( uniquely solvable when a(A) n a(L 22 )=(p) 
with L 22 chosen appropriately stable 


C. Solvability of the Matching Conditions 
Ideal Trajectories: 

1 w* = Aw* + Bu * 
e y =Cw*= 0 

with 


- S 2 z d + G' z TS,z d = QS 2 + G’ z TS,) z d 

s 2 

S,F = AS, + BS 2 = (A + BG* Z T)S X + BS 2 
CS X =0 

S X F = (A + BG * )5, +BS 2 +T6 
CS X = Q9 


Also 

Am = u — u, = u — S 2 z d 
= u~(S 2 z D + G>.) 

= u — S 2 L d (p D — G z z* 

= G]e + (G; - S 2 L d )( p D + G‘ z z- G\ z. 
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■■ Gj (CAw - Q9H D e) + G' z {TAw +e) + h 
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so we can use adaptive control law (32) with gain adaption 
laws (33). 



D. Proof that State Feedback does not Change the 
Minimum Phase of Transfer Function 


Assume a square transfer function (M=P). 
The transmission zeros of (A,B ,C) are 
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=> transmission zeros of (A + B G, B f C) 
are the same as those of (A, Z?, C) # 


(36) 


IV. Implementation of ADTC Theory with Partial 
State Estimation and State Feedback 

A. Simulation Setup 

We implemented the ADTC with partial state feedback to 
the NREL’s 5 MW offshore wind turbine model [7]. This is 
three b laded horizontal axis upwind wind turbine with 63 
meter of rotor radius and 5 MW of rated power. The cut-in, 
rated and cut-out wind speed are 3 m/s, 11.4 m/s and 25 m/s 
respectively. Also, cut-in and rated rotor speeds are 6.9 m/s 
and 12.1 m/s respectively. It has rated Tip Speed of 80 m/s 
with a Tip Speed Ratio of 7.55. 

Depending upon which Degree Of Freedom (DOF) we 
need, it has eleven switches which can be switched on and off 
to add the complexity in the operation. When all DOF’s are 
switched on the wind turbine model consists 31 states. The 
detailed description of this wind turbine can be found in [7]. 

To design the wind speed and partial state estimation of 
wind turbine, the turbine is linearized at constant wind speed 
of 8 m/s with blade pitch held at 0 degree. During the 
linearization Drive Train and Generator DOF switch were 
turned on which gave four states with two states due to the 
Generator DOF and two states due to Drive Train DOF. The 
first state (generator azimuth position) was removed to get 
the three state model of wind turbine since the linearized 
model becomes unobservable when included all four states. 
This three state model is then augmented with wind 
disturbance model given by equation (7) and state estimator 
was designed using equation (9). 

B. Simulation Results 

The ADTC with state estimator and state feedback 
designed is then implemented and simulated in NREL’s 
5MW onshore wind turbine in MATLAB/Simulink platform. 
We also simulated existing fixed gain PID controller [7] to 
compare the performance of ADTC. 


The value of adaptive gains used were yD= 15000 and 
ye=0.0001. We obtained these values of adaptive gain by trial 
and error. 



Time (sec) 

Figure 2. Actual vs estimated wind speed 

Comparison of estimated wind speed with actual step 
wind profile is in figure 3. The wind speed estimator 
combined with the partial state estimator has estimated the 
wind speed to a great extent. 
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Figure 3. Tip Speed Ratio comparison 


The TSR is compared between proposed ADTC and the 
existing fixed gain controller in figure 4. In this simulation, 
the desired TSR is 7.55. The ADTC has better and consistent 
TSR tracking compared to the existing fixed gain controller. 



Figure 4. Comparison of rotor speed (yp) 
and Q* uD 


The components of estimated output tracking error given 
by equation (12) are shown in figure 5. In this simulation the 
value of tracking ratio (Q) used was 1.12. It can be seen that 
the adaptive controller with partial state feedback has done 
excellent job in disturbance tracking. 
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Figure 5. Estimated tracking error 

The Estimated Tracking Error given by equation (12) is 
shown in figure 6. The error is almost zero in steady state. 



The generator speed profile for the adaptive controller and 
the existing fixed controller is in figure 7. In lower wind 
speed region there is discrepancy between the generator 
speed, this is because the existing fixed gain controller the 
generator torque is directly the function of generator speed, 
given by K*Q 2 , where K is some constant, which depends on 
the desired TSR, and Q is the generator speed. But, in 
adaptive disturbance tracking control, there is no direct 
relation between generator speed and torque, instead 
generator torque is manipulated in such a way that the 
generator speed track some percentage of wind speed. 

The generator torque command input produced by both the 
adaptive controller and the fixed controller is shown in figure 
8. The adaptive controller has aggressive torque compared to 
the fixed gain controller. The reason behind this 
aggressiveness is tendency of adaptive controller to track the 
step wind speed. This aggressive behavior can be reduced to 
a great extent by using a low pass filter after estimated wind 
speed. This makes the wind speed smooth rather than step 
resulting the adaptive controller to follow the smooth wind 
speed. 



Figure 7. Generator Torque comparison 
V. Conclusion 

The theory of ADTC with partial state estimator and state 
feedback is proposed to track the TSR in Region II operation 
of large HAWTs, and its stability is also analyzed. The 
estimator/controller is designed based on proposed theory 
and simulated in 5 MW onshore wind turbine. The main 
advantage of this theory is simplicity in design and omission 
of requirement of detailed model of the wind turbine. Also, 
parameters of wind turbines are poorly known which 
prevents design of controller using well established linear 
control theory but adaptive controller performs better in such 
situation. 
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